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ABSTRACT 

Background Mitochondrial DNA (mtDNA) diseases are 
rare disorders whose prevalence is estimated around 1 in 
5000. Patients are usually tested only for deletions and 
for common mutations of mtDNA which account for 
5-40% of cases, depending on the study. However, the 
prevalence of rare mtDNA mutations is not known. 
Methods We analysed the whole mtDNA in a cohort 
of 743 patients suspected of manifesting a mitochondrial 
disease, after excluding deletions and common 
mutations. Both heteroplasmic and homoplasmic variants 
were identified using two complementary strategies 
(Surveyor and MitoChip). Multiple correspondence 
analyses followed by hierarchical ascendant cluster 
process were used to explore relationships between 
clinical spectrum, age at onset and localisation of 
mutations. 

Results 7.4% of deleterious mutations and 22.4% of 
novel putative mutations were identified. Pathogenic 
heteroplasmic mutations were more frequent than 
homoplasmic mutations (4.6% vs 2.8%). Patients 
carrying deleterious mutations showed symptoms before 
16 years of age in 67% of cases. Early onset disease 
(<1 year) was significantly associated with mutations in 
protein coding genes (mainly in complex I) while late 
onset disorders (>1 6 years) were associated with 
mutations in tRNA genes. MTND5 and MTND6 genes 
were identified as 'hotspots' of mutations, with Leigh 
syndrome accounting for the large majority of associated 
phenotypes. 

Conclusions Rare mitochondrial DNA mutations 
probably account for more than 7.4% of patients with 
respiratory chain deficiency. This study shows that a 
comprehensive analysis of mtDNA is essential, and 
should include young children, for an accurate diagnosis 
that is now accessible with the development of next 
generation sequencing technology. 



INTRODUCTION 

Since the first identification of mitochondrial DNA 
(mtDNA) defects responsible for human diseases in 



1988, it has been increasingly clear that mutations 
in the mtDNA represent an important cause of 
neuromuscular disorders/ ^ Nevertheless, today, 
the prevalence of mtDNA disease is still difficult to 
assess for two main reasons.^ First, large series are 
extremely rare because of clinical heterogeneity and 
diagnosis complexity.^ Second, molecular screening 
is usually restricted to detection of mtDNA dele- 
tions and a few^ common mutations w^ithout having 
accurate data on the prevalence of rare mutations. 

Improvements in diagnosis methods led us to 
study a large cohort of patients from the French 
Mitochondrial Disease Network, including both 
children and adults suspected of manifesting a 
mitochondrial disorder. In our cohort, deletions 
and common mutations of mtDNA had already 
been excluded and we performed a comprehensive 
analysis of the mitochondrial genome in order to 
answ^er a number of key questions. 

► What is the prevalence of rare mtDNA muta- 
tions in mitochondrial disorders? 

► What are the phenotypes associated w^ith these 
mutations? 

► Is a comprehensive screening of mtDNA essen- 
tial for optimal diagnostic approach? 

PATIENTS AND METHODS 
Study population 

Over a 3 year period, 743 individuals suspected 
of having a mitochondrial disorder w^ere included 
in the study. All were studied in French referral 
centres by hospital specialists. A standardised chart 
including clinical symptoms, imaging and extensive 
laboratory w^ork w^as completed for all patients (ie, 
familial and clinical history, clinical presentation, 
brain MRI, metabolic screening, mitochondrial 
enzymes studies, histological and molecular ana- 
lyses). Age of onset of clinical symptoms ranged 
from the neonatal period to 74 years of age. The 
patient population w^as based on the follow^ing 
inclusion criteria: (1) cHnical features suggesting 
respiratory chain (RC) dysfunction and/or 
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(2) metabolic screening suggesting RC dysfunction and/or 

(3) isolated or multiple RC complex deficiency and/or (4) histo- 
chemical evidence of mitochondrial abnormality in the 
muscle biopsy and (5) exclusion of mtDNA deletion (s) and 
m.3243A>G, m.8344A>G and m.8993T>C/G point muta- 
tions. Affected individuals presenting w^ith mtDNA depletion or 
clear mendelian family history w^ere also excluded. Among the 
743 patients, 73 satisfied one, 301 two, 275 three and 94 four 
inclusion criteria. Patients w^ere subdivided into three groups 
according to their age at presentation: < 12 months of age, 
1-16 years of age and >16 years of age. 

Blood and tissue samples were obtained after parents of 
affected children and adult patients had given informed consent. 
Authorisations for individual data protection w^ere obtained. 

OXPHOS spectrophotometric measurements 

Spectrophotometric studies of oxidative phosphorylation 
(OXPHOS) complexes and citrate synthase w^ere performed on 
tissue homogenates and fibroblasts, as previously described.^ 

Surveyor and MitoChip analysis 

The entire human mtDNA was first analysed by the Surveyor 
method, as previously described.^ ^ Then, screening w^as per- 
formed w^ith a resequencing chip (MitoChip, GeneChip 
Mitochondrial Resequencing Array 2.0) according to the manu- 
facturer's instructions (Affymetrix, Inc).^ A custom bioinformat- 
ics pipeline was then used for MitoChip analysis, as previously 
described.^ Sequences w^ere compared w^ith the human mtDNA 
consensus sequence, Genbank No J01415.2.^^ 

Variants reported in the MITOMAP database (http://vww. 
mitomap.org/MITOA/LAP) as polymorphisms and present four 
times or more in the mtDB database (http://w^vw.mtdb.igp.uu. 
se/) were considered as non-deleterious and excluded from 
the study. 

Statistical analysis 

Variables 

The variables used refer to clinical data of patients, age at onset 
(<lyear, 1-16 years, > 16 years) and mutation localisation 
(tRNA genes vs protein coding genes). 

Descriptive statistics 

Descriptive statistics w^ere used to characterise the population. 
The unit of analysis w^as the patient and the analyses w^ere done 
on the entire population. All variables w^ere categorical. 

Multiple correspondence analysis and hierarchical ascendant 
classification 

Analyses were performed on the patient population carrying 
pathogenic mutations. The unit of analysis w^as the patient. To 
explore the genotype-phenotype relationships, w^e first used 
multiple correspondence analysis (MCA)^^ and next considered 
the coordinates of the observations on the retained factorial 
axes as new variables used for hierarchical ascendant classifica- 
tion (HAC).^^ MCA allow^s a graphical representation of the 
strength of the association betw^een categories of selected vari- 
ables. The points corresponding to the categories of each vari- 
able (eg, age of onset of disease in three categories) are plotted 
at distances from the origin that are inversely related to the 
number of patients in a given category. Benzecri's method w^as 
used to select the number of axes considered in the analyses. 
MCA w^as built using age at onset and mutation localisation as 
supplementary (non-active) variables as they w^ere the variables 
to explain. The variables w^ith a low frequency (<15%) were 



excluded from the analysis and Cramer's V test was performed 
to identify the variables strongly associated w^ith others, w^hich 
w^ere analysed as non-active variables. To build homogeneous 
clusters of patients, w^e performed HAC based on the Ward 
method. The clustering process w^as plotted as a dendrogram, 
M^ith horizontal branches representing the combination of 
tw^o clusters and vertical branches the degree of dissimilarity 
between combined clusters. Long distances of the vertical seg- 
ments indicated large differences between the clusters. 
Associations betw^een profiles of patients and variables used in 
the model w^ere tested w^ith Fisher's exact test. Test values w^ere 
calculated to measure the association betw^een supplementary 
variables and profiles of patients. We did not use other cluster- 
ing methods, such as K means or latent class analysis, because of 
the small size of the population. We chose hierarchical clustering 
because of the small datasets and the absence of previous 
assumptions about the number of clusters. All statistical analyses 
w^ere performed w^ith SAS V9.3 software (SAS Institute, Gary, 
North Carolina, USA). A two tailed p<0.05 w^as considered stat- 
istically significant. 

RESULTS 

The 743 patients were subdivided into 395 men and 348 
w^omen (ratio 1.13 : 1) w^ith an average age of 25 years at the 
time of the study. Clinical data are shown in the online 
supplementary table SI for the entire population, and depend- 
ing on age of onset. Early clinical symptoms were seen during 
the first year of Hfe for 356 (48%) patients, between 1 year and 
16 years for 192 (25.8%) patients and in adulthood for 
the remaining 180 (24.2%) patients. In 15 (2%) of 743 cases, 
age of onset was unknown. Neurological presentations were the 
most frequent (72.7% of total population). Before 1 year, the 
most common symptom was psychomotor retardation (55.6%) 
whereas cerebellar ataxia was commonly seen after 1 year of age 
or during adulthood (37% and 25%, respectively). Epilepsy and 
myoclonies were also very common (32.2%) whatever the age 
of presentation. Muscular involvement was among the second 
rank of clinical disorders, with chronic progressive external 
opthalmoplegia (CPEO) mainly seen during childhood and 
adulthood (27.1% and 32.2%, respectively). Postnatal growth 
failure was highly represented before 1 year (29.5%). 
Hypertrophic cardiomyopathy, and liver and kidney deficiencies, 
were also mainly associated with onset before 1 year (11%, 
27.8% and 11.5%, respectively). Optic atrophy was most fre- 
quent when the disorder started during childhood (11.5%) 
whereas peripheral neuropathy was present in 24.4% of patients 
with adult onset. Paraclinical investigations are presented in the 
online supplementary table S2. Brain MRI was abnormal in 
82.1% of cases before 1 year of age, with a diagnosis of Leigh 
or Leigh-like syndrome in 43.1% of patients. In patients with 
adult onset, white matter involvement was the most common 
pattern (43.8%) whereas posterior fossa atrophy/hypoplasia 
was the most frequent when onset was during childhood 
(29.5%). A lactate peak was detected in 60.6% of cases by 
MRI spectroscopy, mainly in onset before 1 year (68.1%). 
Hyperlactacidaemia was present in about 50% of patients but 
hyperlactatorachia was the most common metabolic abnormality 
(75.3%), mainly during childhood (83.3%) (see online supple- 
mentary table S2). Results obtained from tissue biopsies, mainly 
skeletal tissue, are presented in the online supplementary 
table S3. The percentage of ragged red fibres (RRF) and COX 
negative fibres increased with age of onset, reaching 48.9% and 
56.9%, respectively, after 16 years of age. In muscle, the bio- 
chemical activity of RC complexes was significantly reduced in 
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59.4% of total biopsies. In liver tissue, a RC deficiency was seen 
in 53.9% of cases before 1 year of age. 

The presence of deletions and common point mutations 
(m.3243A>G, m.8344A>G and m.8993T>C/G) of mtDNA 
had previously been ruled out. The first part of the study aimed 
to search for heteroplasmic mutations by Surveyor analysis. 
After eliminating the know^n polymorphisms (list available on 
request) and, among new^ sequence variations, the 52 synonym- 
ous variants and the 28 nucleotide substitutions in the D loop 
region (list available on request), a total of 82 different hetero- 
plasmic sequence variations corresponding to putative or know^n 
pathogenic mutations w^ere retained (see online supplementary 
table S4)}^~^^ The 56 different variants corresponding to puta- 
tive mutations w^ere not reported in the MITOAdAP database 
and were present less than four times in the mtDB database. 
Two variants only were found tw^ice. The m.l938A>T variant 
in MTRNR2 was found in a child who presented w^ith psycho- 
motor delay, epilepsy, ptosis, deafness and w^hite matter lesions 
on brain MRI. The second patient w^as an adult presenting 
M^ith epilepsy, cerebellar ataxia and stroke-like episodes. 
Neuroimaging revealed patterns of Leigh/MELAS (mitochon- 
drial encephalomyopathy, lactic acidosis and stroke-like epi- 
sodes) overlap syndrome. In both cases, biochemical analysis 
revealed no RC deficiency on muscle biopsy. The second 
variant, m.l3604G>C, w^as in the MTDNS gene. The first child 
died at birth and presented a complex I (CI) deficiency in the 
liver. The second child presented w^ith myopathy, CPEO and 
deafness. Muscle biopsy show^ed RRF and COX- fibres but bio- 
chemical analysis was not done. 

Tw^enty-six different sequence variations w^ere already know^n 
as pathogenic mutations. The characteristics of the 34 corre- 
sponding patients (34/743, 4.6%) are presented in table 1. Four 
mutations (m.9185T>C, m.l0191T>C, m.l2706T>C and 
m.l3514A>G) were found in two patients, m.l4487T>C in 
three patients and m.l3513G>A in four patients. Recurrent 
mutations w^ere responsible for Leigh syndrome in these 
patients, except for one, and encoded subunits of CI, except the 
m.9185T>C variant located in the MTATP6 gene (table 1). This 
w^ork expands the clinical spectrum of mtDNA disorders, w^ith 
the example of m.636A>G and m.l2236G>A mutations that 
had been previously reported in non-syndromic hearing loss and 
that are associated w^ith progressive multisystemic disorders in 
our cohort. 

In 74% of cases (25 of 34patients), symptoms started before 
16 years of age (figure 1). When onset w^as before 1 year, all 
patients presented a multisystemic disorder w^ith abnormal brain 
MRI, hyperlactataemia and elevation of the lactate to pyruvate 
ratio. Leigh syndrome w^as found in 85% of cases (six of seven 
patients), and usually associated w^ith a CI deficiency in muscle 
(four of six patients w^ith muscle biopsy, 66.7%). When the 
onset was between 1 and 16 years of age, ocular symptoms were 
present in 44.4% of patients (8/18). In adult onset, cerebellar 
ataxia w^as found in 44.4% of cases (four of nine patients) and 
biochemical analysis identified a multiple RC deficiency in 
muscle in 42.9% (three of seven patients w^ith muscle biopsy). 
In agreement w^ith the mitochondrial enzyme data, mutations 
were mainly located in CI genes before 1 year (five of seven 
patients, 71.4%) and in tRNA genes in adult onset (six of nine 
patients, 66.7%) (figure 1). Statistical analysis confirmed these 
results. MCA, used as a step of data reduction, retained 12 
active variables: affected maternal relatives, psychomotor retard- 
ation, regression, hypotonia, stroke-like episodes, epilepsy, 
movement disorders, CPEO, optic atrophy, deafness, postnatal 
groM^th retardation and Leigh syndrome. The follow^ing variables 



were strongly correlated w^ith others and w^ere added as non- 
active: cerebellar ataxia, w^hite matter involvement and posterior 
fossa atrophy or hypoplasia. HAC results show^ed three distinct 
profiles (figure 1). The first profile of patients w^as associated 
w^ith early onset disease (<lyear) and mutations in protein 
coding genes. The follow^ing variables w^ere overrepresented in 
this group of patients (p<0.05): psychomotor retardation, 
hypotonia, optic atrophy, postnatal grow^th retardation and 
Leigh syndrome. The second profile of patients was associated 
M^ith late onset (> 16 years) and mutations in tRNA genes. In 
this profile, affected maternal relatives were common (36%) 
(p = 0.13). We did not find any association for the last profile 
of patients w^ith a disease onset during childhood (1-16 years). 
A dendrogram for cluster model is illustrated in the online 
supplementary figure SI. 

The second part of the study was to identify homoplasmic 
mutations using MitoChip analysis in patients w^ithout patho- 
genic heteroplasmic variants. Know^n polymorphisms and 180 
novel variants (140 synonymous and 40 in the D loop) were 
excluded (list available on request). The online supplementary 
table S5 show^s the 120 different homoplasmic variants found in 
131 patients, w^hich correspond to potential and deleterious 
mutations. Among the 104 variants corresponding to putative 
pathogenic mutations, four w^ere present more than once (see 
online supplementary table S5). The m.3228_3229 insA variant 
in MTRN2 was found in two adults w^ith a mitochondrial myop- 
athy and in one child w^ith Leigh syndrome. The m.8654T>C 
(p.Ile43Thr) variant in MTATP6 was associated w^ith exercise 
intolerance and multiple RC deficiency in muscle of an adult 
patient. The second affected individual w^as a child w^ith Leigh 
syndrome. None show^ed complex V deficiency. The last two 
variants w^ere in MTCYB. Two adult patients, presenting w^ith 
exercise intolerance, carried the m.l5122A>G (p.Thrl26Ala) 
variant w^ith classical mitochondrial features on muscle biopsy 
and CIII deficiency in only one of two. The m.l5774T>C 
(p.Val343Ala) variant was associated w^ith tw^o early onset pre- 
sentations (<1 year). The first child presented w^ith Toni Debre 
Eanconi syndrome and the second one had a neurological pres- 
entation (psychomotor delay, hypotonia and epilepsy). In both 
cases, RRE fibres were observed on muscle biopsy but we found 
no RC deficiency. 

Sixteen different sequence variations were already know^n 
as pathogenic mutations and the characteristics of the 21 corre- 
sponding patients (2.8%) are presented in table 2. The 
m.l555A>G mutation, only seen in a homoplasmic state, w^as 
present in three patients w^ith a multisystemic disorder; only one 
was deaf. Two other mutations were present more than once, 
and both had been previously found in our study in the hetero- 
plasmic state (m.4317delA, m.l3513G>A). 

Ocular involvement w^as observed in tw^o-thirds of patients 
and deafness in half of affected individuals during childhood 
(figure 2). HoM^ever, the population size w^as too small to 
perform statistical studies. 

Analysis of all (hetero- and homoplasmic) deleterious muta- 
tions confirmed the role of rare mtDNA mutations in paediatric 
presentations (onset < 16 years for 37 of 55 patients, 67.3%) 
(figure 3). Neuromuscular symptoms were clearly predominant 
in the 55 patients carrying pathogenic mutations. Liver failure 
was associated w^ith onset before 1 year, ocular involvement was 
present in 50% of patients (12/24) in childhood presentations 
and diabetes mellitus and deafness w^ere found in 20.8% of 
cases in adult onset disorders (5 of 24 patients). Statistical 
analysis confirmed previous findings. The same active and non- 
active variables were retained by MCA. Results of HAC on the 
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Table 1 Characteristics of patients carrying pathogenic mutations identified by Surveyor analysis 
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55 patients showed three distinct profiles (figure 3). Early onset 
disease and mutations in protein coding genes were associated 
with a profile of patients characterised by psychomotor retard- 
ation, regression, stroke-like episodes, movement disorders, 
Leigh syndrome and white matter involvement. Late onset was 
associated with mutations in the tRNA genes. Affected maternal 
relatives were also found predominantly in this second profile 
(33%), despite being non-significantly associated (p = 0.18). 
Again, we found no pattern associated with an onset during 
childhood. A dendrogram for cluster model is illustrated in the 
online supplementary figure S2. 



DISCUSSION 

In OXPHOS disorders, identification of a causative mutation(s) 
is critical to confirm the diagnosis and to propose accurate 
genetic counselling and prenatal diagnosis, but molecular ana- 
lysis is problematic due to the large genetic heterogeneity of 
mtDNA related disorders. Primary mtDNA defects include large 
rearrangements and point mutations. In clinical practice, most 
suspected patients do not have their entire mitochondrial 
genome sequenced but are usually tested only for deletions and 
a small number of common mutations. Several groups have 
reported that testing for common mtDNA defects in adults with 
suspected mitochondrial disease can have quite a high diagnostic 
yield, ranging from 10% to 40% of patients tested.^ ^ Testing 
for common mtDNA mutations has a much lower diagnostic 
yield in children, ranging from 2-4% to 25%.^^ Recently, 
looking for 10 mtDNA point mutations in 3168 neonatal cord 
blood samples revealed that at least 1 in 200 healthy humans 
harbours a pathogenic mtDNA mutation, although at low 
mutant heteroplasmy.^^ However, to date, the rate of rare muta- 
tions is totally unknown. To answer this question, we analysed 
the whole mitochondrial genome of 743 patients suspected of 
manifesting a mitochondrial disease after excluding deletions 
and the most common mutations. 

A total of 55 patients (7.4%) harboured pathogenic mutations 
responsible for an onset disorder before adulthood in 67% of 
cases and before 1 year of age in 24% of cases. Previous studies 
have suggested that 50% or more of all patients have adult 
onset. However, in these studies, the most common mutations 
(mainly deletions and m.3243A>G) that were mainly found 
in adult presentations (CPEO, maternally inherited deafness 
and diabetes (MIDD), etc) were not excluded. The French 
Mitochondrial Disease Network includes approximately 90% of 
patients investigated for mitochondrial disease in France. The 
743 patients included in the study were doubly selected because 
they were highly 'suspect' of RC disorder, and they did not 
carry common mtDNA mutations. Our results suggest that in 
the French population, the clinical presentations that are not 
associated with common mutations begin mainly before adult- 
hood. Further studies using next generation sequencing technol- 
ogy will be necessary to confirm this point. 

Assuming that heteroplasmic and homoplasmic mutations 
may have different effects, we also analysed both populations 
separately. As expected, heteroplasmic mutations were more fre- 
quent than homoplasmic mutations (4.6% vs 2.8%). Symptoms 
appeared before the age of 16 years in 75% of patients carrying 
heteroplasmic pathogenic mutations and in 58% of patients 
carrying homoplasmic pathogenic mutations. If a number of 
factors have led to an exaggerated perception that the mtDNA 
genome is rarely involved in children affected by mitochondrial 
disorders, our results show that rare mtDNA mutations are a 
non-negligible cause of paediatric RC deficiency because 6.75% 
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Figure 1 Analysis of the population harbouring heteroplasmic pathogenic mitochondrial DNA (mtDNA) mutations according to age at onset, 
n, number of patients. (A) Age at onset (in years). (B) Clinical characteristics. CNS, central nervous system; D, deafness; DM, diabetes mellitus. 
GR, growth retardation (intrauterine or postnatal); HCM, hypertrophic cardiomyopathy; L, liver involvement; Mu, muscle; OC, ocular involvement 
(optic atrophy or pigmentary retinopathy); PN, peripheral neuropathy; Re, renal involvement. (C) Respiratory chain analysis by spectrophotometry. 
CI, CV, respiratory chain deficiency in complex I, in complex V; mul, multiple respiratory chain deficiency; N, normal. (D) Localisation of mtDNA 
mutations. CI, CI 1 1, CV, genes encoding subunits of complexes I, III, V; tRNA, genes encoding tRNA. (E) Statistical analyses with hierarchical 
ascendant classification of pathogenic heteroplasmic mutations found in the 34 patients. The upper right quadrant is complementary to the bottom 
left quadrant. The numbers on the axes are the scores that represent the contribution of each feature to the overall inertia. Our model (axes 1-3) 
explains 48% of total inertia. The circles represent the variables. Their size is proportional to the quality of representation of the variables on the 
plane (axes 1 and 3). Patients with disease onset <1 year are clustered with mutations in protein coding genes, psychomotor retardation, hypotonia, 
optic atrophy, postnatal growth retardation and Leigh syndrome. Patients with disease onset >1 6 years are clustered with mutations in tRNA genes 
and affected maternal relatives. 
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Table 2 Characteristics of patients with pathogenic mutations identified by MitoChip analysis 



Gender/age at OH-but/ 



Gene 
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onset 


Clinical symptoms 


Brain MRI 
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RC 
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ND 
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N 
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LUa— 


M 
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nn.3460A>G 


M>16 


CNS (A), PN, Mu (ptosis), HCM, 


Leigh 


ND 


N 


ND 


RRF, COX-, 


iCI (Mu) 
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m.3890G>A 
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t (B, CSF) 


t 


ND 


Lipidosis 


N (Mu) 








OA 
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m.4300A>G 


F>16 


HCM 


ND 


ND 


ND 


ND 
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ND 
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ND 


ND 
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M 
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111. J J'+UVJ>M 
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IN 


IMU 
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1 ri\/ (M\\\ 
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lipidosis 




MTATP6 


m.9185T>C 


F>16 
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ND 


ND 


ND 


ND 


N 


N (Mu) 




m 1 ^G>A 


1 ^IVI^ 1 u 


CNS (fn F M-IH) PN Mu 

v_iMj \tii, t, IVI in/, r iM, IVIU 


r~prphpllar atrnnhw 
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ND 


Mitnrhnnrlrial 

IVIl LUV_I lUI lUI iOI 


If! (Mu) 
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(SM), HCM, PR, OA, D 


Stroke-like 
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1<M<16 


CNS (PMR, A), Mu (PEO), OA, 


ND 


ND 


ND 


ND 


ND 


ND 








lU-G, PN-G 














MTND6 


m.14459G>A 


F<1 


CNS (R) 


Leigh 


N(B) 


ND 


ND 


ND 


N (Mu), iCI (L), 




















iCI+IV (F) 






M<1 


CNS (PMR, Dys), Mu (ptosis) 


Leigh 


t (B, CSF) 


ND 


ND 


ND 


iCI (Mu, F) 




m.14487T>C 


1<M<16 


CNS (PMR, S-l, PS), Mu (ptosis) 


Leigh 


t (CSF) 


N 


N 


N 


iCI (Mu) 


MTCYB 


m.15234G>A 


1<M<16 


CNS (R, PMR, E, S-l, M-IH), Mu 


Leigh, Stroke-like, 


t (B, CSF) 


t 


t 


RRF 


iCI, Clll, CIV 








(PEO, SM), D, cataract 


leukodystrophy. 










(M), N (F) 



i, decreased; t, increased; B, blood; CI, II, III, IV, V, respiratory chain complexes; CNS, central nervous system (A, cerebellar ataxia; Dys, dystonia; E, epilepsy; En, encephalopathy; 
MD, movement disorders, M-l H, migraine-like headaches; PMR, psychomotor retardation; PS, pyramidal syndrome; R, regression; S-l, stroke-like episode); COX-, COX negative fibres; 
D, deafness; DM, diabetes mellitus; F, fibroblasts; HCM, hypertrophic cardiomyopathy; lU-G, intrauterine growth failure; L, liver involvement; L/P, lactate/pyruvate; Mu, muscle (H, 
hypotonia; PEO, progressive external ophthalmoplegia; SM, skeletal myopathy); N, normal; ND, not determined; OA, optic atrophy; OH-but/Ac-Ac, hydroxy-butyrate/aceto-acetate; 
PN, peripheral neuropathy; PN-G, postnatal growth failure; PR, pigmentary retinopathy; RC, respiratory chain; RRF, ragged red fibres; WM, white matter involvement. 



of patients with onset before adulthood (37 of 548) carried rare 
pathogenic mtDNA mutations. 

One of the main results of this study was that the location of 
mtDNA mutations influenced the age of onset of OXPHOS dis- 
eases. We found a significant correlation between mutations in 
protein coding (mainly CI) genes and early onset disorder 
(<lyear) on the one hand, and between mutations in tRNA 
genes and adult onset on the other hand; the hierarchical clus- 
tering allowing unambiguous class assignments of patients with 
early and late onset disease. These results suggest that mutations 
in tRNA genes are much better tolerated than those that directly 
affect the functioning of CI whose defect is the most frequently 
observed. It was previously thought that a mtDNA origin 
accounted for only 5-10% of paediatric cases. However, a 
recent study of 109 children with isolated CI deficiency identi- 
fied pathogenic mtDNA mutations in 29% of cases.^^ In our 
study, mutations in MTND genes were found in 11 of 75 
patients (14.7%) with an isolated CI deficiency before the age 
of 16 years but this percentage is likely underestimated because 
other deleterious mutations were likely present among the 44 
variants corresponding to putative mutations that we identified 



in MTND genes. Among CI genes, MTNDS and MTND6 corre- 
sponded to 'hotspots' for disease causing mutations, being, 
respectively, found mutated in 11 (20%) and seven (12.7%) of 
55 patients , with Leigh syndrome accounting for the large 
majority of associated phenotypes. 

A total of 167 patients (22.4%) harboured putative patho- 
genic mutations, with homoplasmic novel variants being 
more frequent than heteroplasmic variants (14.6% vs 7.8%). 
Interpretation of the clinical significance of rare variants is chal- 
lenging due to the highly polymorphic nature of mtDNA. The 
existing scoring systems to evaluate the pathogenicity of a 
mtDNA variant consider several parameters, including single 
fibre PCR or cybrid studies, multiple family reports, segregation 
within the family, measurable biochemical defect, evolutionary 
conservation, etc.^^ Collection of data and tools to perform 
functional studies is in progress but we currently do not have 
the information needed to precisely define the status of the 160 
different variants identified. Nevertheless, 32 patients (4.3%) 
carried novel variants in protein coding genes with high or 
medium amino acid conservation and a damaging and not toler- 
ated status according to our in silico analysis using PolyPhen 2 
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Respiratory chain analysis Localisation of mtDNA mutations 

in muscle 

Figure 2 Analysis of the population harbouring homoplasmic pathogenic mitochondrial DNA (mtDNA) mutations according to age at onset, 
n, number of patients. (A) Age at onset (in years). (B) Clinical characteristics. CNS, central nervous system; D, deafness; DM, diabetes mellitus; 
GR, growth retardation (intrauterine or postnatal); HCM, hypertrophic cardiomyopathy; L, liver involvement; Mu, muscle; OC, ocular involvement 
(optic atrophy or pigmentary retinopathy); PN, peripheral neuropathy; Re, renal involvement. (C) Respiratory chain analysis by spectrophotometry. 
CI, cm, CIV, respiratory chain deficiency in complex I, complex III, complex IV; mul, multiple respiratory chain deficiency; N, normal. (D) Localisation 
of mtDNA mutations. CI, CI 1 1, CV, genes encoding subunits of complexes I, III, V; rRNA, genes encoding rRNA; tRNA, genes encoding tRNA. 



and SIFT, respectively (see online supplementary tables S4 
and S5).^^ In 15 of 32 patients , the identified variant w^as corre- 
lated M^ith a corresponding biochemical defect (not show^n). 
Regarding tRNA coding genes, mutations w^ere detected in 
18 of 22 genes and were found spread all over the structural 
domains of the corresponding tRNAs. Seventeen of the 36 
new^ly detected variants affect highly conserved nucleotides 
(table 3). Despite the fact that conservation is often poor in 
picking pathogenicity,^^ this is a criterion that may be taken 
into account to prioritise future functional studies allow^ing 
assessment of tRNA variants using scoring systems. 

In conclusion, w^e have reported the first study to determine 
the prevalence of rare mtDNA mutations in RC disorders. 



Although further analyses are needed to define precisely the 
pathogenic or polymorphic nature of all of the novel variants 
identified, our results show^ that rare mtDNA mutations 
account for more 7.4% of patients w^ith mitochondrial disor- 
ders. Careful selection of the population studied and the 
choice of tissues concerned by molecular analysis, w^hich was 
carried out mainly in muscle (88% of cases), may explain this 
mutation rate, probably higher if one takes into account some 
new^ variants are probably deleterious. In this cohort of 
patients, w^e excluded those carrying a deletion or common 
mtDNA mutations. It is interesting to note that w^ithin the 
French Mitochondrial Disease Netw^ork, we have a rate of 
8% for these common abnormalities of the mitochondrial 
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Figure 3 Analysis of the pooled populations harbouring heteroplasmic or homoplasmic pathogenic mitochondrial DNA (mtDNA) mutations 
according to age at onset, n, number of patients. (A). Age at onset (in years). (B) Clinical characteristics. CNS, central nervous system; D, deafness; 
DM, diabetes mellitus; GR, growth retardation (intrauterine or postnatal); HCM, hypertrophic cardiomyopathy; L, liver involvement; Mu, muscle; OC, 
ocular involvement (optic atrophy or pigmentary retinopathy); PN, peripheral neuropathy; Re, renal involvement. (C) Respiratory chain analysis by 
spectrophotometry. CI, Clll, CIV, CV, respiratory chain deficiency in complex I, complex III, complex IV, complex V; mul, multiple respiratory chain 
deficiency; N, normal. (D) Localisation of mtDNA mutations. CI, Clll, CV, genes encoding subunits of complexes I, III, V; rRNA, genes encoding rRNA; 
tRNA, genes encoding tRNA. (E). Statistical analyses with hierarchical ascendant classification of pathogenic mutations found in the 55 patients. Our 
model (axes 1-3) explains 45% of total inertia. The circles represent the variables. Their size is proportional to the quality of representation of the 
variables on the plane (axes 1 and 2). Patients with disease onset <1 year are clustered with protein coding gene mutations, psychomotor 
retardation, stroke-like episodes, regression, movement disorders, Leigh syndrome and white matter involvement. Patients with late onset are 
clustered with tRNA gene mutations and affected maternal relatives. 
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Table 3 Novel variants in tRNA genes and tentative prediction of impact according to the level of nucleotide conservation of the affected 
position 





Sequence 
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Structural domains in 


Position in 


Change in 


Nucleotide 




Gene 


variation 


reading 


tRNA 


tRNA 


tRNA 


conservation 


Tentative prediction 


MTTF 


fil.D/y 1 >L 


DirGCt 


Acc stem 


0 

J 


1 >L 


yU /0<j{<iiUU /o, 1 or n 


rossiuly patnoloQical 




rn.jyzL> i 
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D loop 


lb 


L> 1 


<jU 70 


Polymorphism 




iTi. jyj 1 >L 


Direct 


D loop 


1 7 


1 >L 


<jU /o 


Polymorphism 




m.645A>G 


Direct 


Acc stem 


71 


A>G 


90%<x<100%, allR 


Possibly pathological 


MTTV 


m.1628C>T 


Direct 


Anticd stem 


29 


C>T 


<50% 


Polymorphism 




m.1645A>G 


Direct 


Variable region 


46 


A>G 


100% A 


Possibly pathological 


MULI 


m.3258T>C 


Direct 


Anticd stem 


27 


T>C 


90%<x<100%, allR 


Possibly pathological 


MUQ 


m.4354C>T 


Inverse 


Variable region 


48 


G>A 


<50% 


Polymorphism 




m.4394C>A 


Inverse 


Acc stem 


7 


G>T 


<50% 


Polymorphism 


MUM 


m.4449G>A 


Direct 


7"5feA77 


52 


G>A 


100%G 


Possibly pathological 


MTTW 


m.5527A>G 


Direct 


D loop 


16 


A>G 


50%<x<90%, Y or R 


Probable 
polymorphism 




m.5542C>T 


Direct 


Anticd loop 


32 


C>T 


100%C 


Possibly pathological 




m.5560G>A 


Direct 


T stem 


51 


G>A 


50%<x<90%, all R 


Probable 
polymorphism 




m.5574T>C 


Direct 


Acc stem 


68 


T>C 


50%<x<90%, all Y 


Probable 
polymorphism 


Mm 


m.5605A>G 


Inverse 


T loop 


55 


T>C 


90%<x<100%, YorR 


Possibly pathological 




m.5608C>T 


Inverse 


T stem 


52 


G>A 


50%<x<90%, all R 


Probable 
polymorphism 


MTTN 


m.5663C>T 


Inverse 


Acc stem 


67 


G>A 


50%<x<90%, all R 


Probable 
polymorphism 




m.5690A>G 


Inverse 


Anticd stem 


40 


T>C 


100%T 


Possibly pathological 


MTTC 


m.5793A>G 


Inverse 


T loop 


60 


T>C 


<50% 


Polymorphism 




m.5800A>G 


Inverse 


Anticd loop 


32 


T>C 


90%<x<100%, ally 


Possibly pathological 


MTTY 


m.5865T>C 


Inverse 


Anticd stem 


31 


A>G 


90%<x<100%, allR 


Possibly pathological 


MUSI 


m.7455A>G 


Inverse 


T stem 


64 


T>C 


90%<x<100%, ally 


Possibly pathological 


MUD 


m.7567C>T 


Direct 


T loop 


54 


C>T 


50%<x<90%, Y or R 


Probable 
polymorphism 


MTTK 


m.8306T>C 


Direct 


D-stem 


12 


T>C 


100%T 


Possibly pathological 


MUH 


m.12168G>A 


Direct 


Anticd loop 


34 


G>A 


100%G 


Possibly pathological* 




m.12173T>C 


Direct 


Anticd stem 


39 


T>C 


50%<x<90%, Y or R 


Probable 
polymorphism 




m.12174C>T 


Direct 


Anticd stem 


40 


C>T 


90%<x<100%, ally 


Possibly pathological 




m.12197_12198lnsC 


Direct 


T stem 


64-65 


Ins C 


1 


Probable 
polymorphism 


MTfLZ 


m.12324C>T 


Direct 


7"5feA77 


61 


C>T 


100%C 


Possibly pathological 


Mm 


m.14697C>T 


Inverse 


T stem 


50 


G>A 


50%<x<90%, Y or R 


Probable 
polymorphism 


Mm 


m.15916T>C 


Direct 


Anticd stem 


31 


T>C 


50%<x<90%, all Y 


Probable 
polymorphism 




m.i5926C>T 


Direct 


Anticd stem 


41 


T>C 


90%<x<i00%, all Y 


Possibly pathological 




m.15936A>T 


Direct 


T stem 


52 


A>T 


50%<x<90%, Y or R 


Probable 
polymorphism 




m.15947A>G 


Direct 


Acc stem 


67 


A>G 


50%<x<90%, Y or R 


Probable 
polymorphism 


Mm 


m.15977C>T 


Inverse 


T stem 


50 


G>A 


50%<x<90%, all R 


Probable 
polymorphism 




m.15992A>T 


Inverse 


Anticd loop 


34 


T>A 


100%T 


Possibly pathological* 



Nucleotide conservations are retrieved from the compilation of tRNA genes of each specificity taking into account 150 mammalian mitochondrial genomes (http://mamit-trna.u-strasbg. 
fr/).^^ Four levels of conservation are considered. Whenever a nucleotide is conserved in more than 90% of the sequences, the tentative prediction is that the probability of having a 
structural or functional role is high. Its mutation may thus possibly have pathological consequences. Possible pathological variants are shown in italic. 

*Position 34 corresponds to the first nucleotide of the tRNA anticodon triplet. The mutations concerning both MTTH and Mm do not affect the codon reading. However, they may 
affect the efficiency of aminoacylation or of post-transcriptional modification, leading to less efficient translation. 
Acc stem, acceptor stem; Anticd loop, anticodon loop; Anticd stem, anticodon stem; R, purine base; Y, pyrimidic base. 



genome. The sum of the rates obtained in our netw^ork for 
both common and rare mutations show^s how^ the complete 
screening of mtDNA is essential. It should be performed for 
an optimal diagnosis, including in young children, and w^ill 



now be accessible w^ith the development of next generation 
sequencing technology that provides a one step approach in 
detecting common and uncommon point mutations, as w^ell as 
deletions. 
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